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An accurate knowledge of the dark matter distribution in the Milky Way is of crucial importance 
for galaxy formation studies and current searches for particle dark matter. In this paper we set new 
dynamical constraints on the Galactic dark matter profile by comparing the observed rotation curve, 
updated with a comprehensive compilation of kinematic tracers, with that inferred from a wide range 
of observation-based morphologies of the bulge, disc and gas. The generalised Navarro-Frenk-White 
(NEW) and Einasto dark matter profiles are fitted to the data in order to determine the favoured 
ranges of local density, slope and scale radius. For a representative baryonic model, a typical local 
circular velocity vq = 230km/s and a distance of the Sun to the Galactic centre Ro = 8kpc, we 
find a local dark matter density po — 0.420i/Q Q^g {2a) ± 0.025 GeV/cm^ (po = 0-420 /Iq q 2 i ± 
0.026 GeV/cm^) for NEW (Einasto), where the second error is an estimate of the systematic due 
to baryonic modelling. Apart from the Galactic parameters, the main sources of uncertainty inside 
and outside the solar circle are baryonic modelling and rotation curve measurements, respectively. 
Upcoming astronomical observations are expected to reduce all these uncertainties substantially 
over the coming years. 


I. INTRODUCTION 


The observation of flat rotation curves in spiral galax¬ 
ies was historically one of the cornerstones upon which 
the case for dark matter was built (e.g. m-. for more 
recent studies, see e.g. [SHS])- In our Galaxy, a spiral 
itself, measuring the rotation curve remains however a 
daunting task due to our position inside the stellar disc. 
As a consequence, the knowledge of the Galactic distri¬ 
bution of dark matter is affected by sizeable uncertain¬ 
ties, especially in the inner regions, where baryons dom¬ 
inate the gravitational potential. This unfortunate cir¬ 
cumstance hinders the study of the formation history of 
the Milky Way and the interpretation of particle dark 
matter searches (see e.g. [MIS]). 

Several techniques have however been designed to pro¬ 
vide gravitational constraints on the dark matter distri¬ 
bution in the Galaxy, and in particular in the solar neigh¬ 
bourhood. They fall into two classes: local methods m- 
[2U] (see also Enna), based on a Jeans-Poisson analysis 
of the kinematics of local stars, and global methods [231- 
[52] , which involve the mass modelling of the Galaxy as a 
whole and use global dynamical constraints. Each class of 
methods brings complementary information on the dark 
matter distribution and both suffer from systematics. For 
example, the precise measurements sometimes provided 
by global methods come at the expense of enforcing a 
preassigned shape to the dark matter prohle in the inner 
Galaxy, where the effect of baryons is potentially large. 
This in turn leads to sizeable systematics in the determi¬ 


nation of the dark matter profile due to the uncertainties 
on the distribution of baryons. 

In Ref. [25], we have applied a global method to show 
that interesting constraints on the dark matter distribu¬ 
tion in our Galaxy can be achieved by comparing the 
rotation curve inferred from baryonic models with actual 
rotation curve measurements. In light of our recent work 
[33] . here we set out to revisit the previous analysis in 
Ref. [29] and improve it with three key ingredients: 


i) a new and more complete compilation of rotation 
curve data, including 2780 measurements from gas 
kinematics, star kinematics and masers (Sec. 0; 


ii) a wide range of baryonic models, including virtually 
all data-based morphologies for the bulge, disc and 
gas (Sec. |n]); and 


iii) a procedure based on a tw o-dimensional applied 
to unbinned data (Sec. III). 


The results, which supersede tho se of our previous 
work [29] . are presented in Sec. IV in the form of con¬ 
straints on the dark matter profile. Our conclusions are 
given in Sec. jV] 


II. OBSERVATIONS 

Our analysis of the dark matter distribution in the in¬ 
ner Galaxy relies upon two observational inputs: the ro¬ 
tation curve and the spatial distribution of stars and gas. 
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Let us begin with the former. Rotation curve measure¬ 
ments in the Milky Way have been available for several 
decades now, but the corresponding literature is rather 
scattered. In Ref. [33] we have assembled a new, com¬ 
prehensive compilation of kinematic tracers across the 
Galaxy that comprises: 

• gas kinematics, including HI terminal velocities 

[34H36] . HI thickness [37], CO terminal velocities 
[35MT] , HII regions [M] and giant molecu¬ 

lar clouds [35] : 

• star kinematics, including open clusters [46], plan¬ 
etary nebulae m. classical cepheids [351 IlH] and 
carbon stars [50] [51]; and 

• masers [53l456| . 

This literature survey is particularly exhaustive in the 
range of Galoctocentric radii i? = 3 — 20 kpc. It inten¬ 
tionally leaves out tracers with relevant random motions, 
asymmetric drift or kinematic distances only in an effort 
to track Galactic rotation as reliably as possible. From 
the references listed above we additionally exclude indi¬ 
vidual objects close to the Galactic centre or anti centre, 
with incomplete data or in some way signalled as sus¬ 
picious. The final compilation consists of 2780 objects 
across the range R = 0.5 — 25 kpc, of which 2174, 506 
and 100 are contributed by gas kinematics, star kine¬ 
matics and masers, respectively. Each object is char¬ 
acterised by its Galactic coordinates {i, b), heliocentric 
distance d, Galactocentric radius R = (c?^ cos^ b + R^ — 
2Rodcosb cos£)^/^ (with Rq the distance of the Sun to 
the Galactic centre) and line-of-sight velocity in the lo¬ 
cal standard of rest (LSR) u}°®. Uncertainties on d and 
are carefully extracted or estimated from the orig¬ 
inal references, whereas those on £ and b are neglected 
since they are by far sub-dominant in all instances. The 
angular velocity uJc{R) and corresponding uncertainties 
are then found by inverting (and propagating) the well- 
known expression 

= (-^oWc - Vo) cos bsini, (1) 

where vq = Vc{Rq) is the local circular velocity. Note 
that we use Wc instead of Vc = RoJc to avoid the cor¬ 
relation between the errors of R and Vc- In the end, 
given Ro, vg and a local standard of rest, each object 
provides an independent constraint on R and uJc[R). 
Throughout the work, unless otherwise specified, we use 
Ro = 8 kpc, Vo = 230km/s and the peculiar solar motion 
{U,V,W)q = (11.10,12.24, 7.25) km/s [57]. The release 
of a user-friendly tool to retrieve all kinematic data used 
here is planned for the near future [55] . 

As for baryons, despite the improvements in the 
amount and quality of observations, their distribution in 
the Milky Way is still not known to high accuracy. In 
order to address the uncertainty on the density profile of 
stars and gas, we have carried out in Ref. |33j an exten¬ 
sive literature survey (cf. Tab. |l]) to bracket the allowed 



model 

specification 

Ref. 


1* 

exponential E2 



2 

gaussian G2 

[591 


3 

gaussian plus nucleus 

m 

bulge 

4 

truncated power law 

m 


5 

power law plus long bar 

m 


6 

truncated power law 

EU 


7 

double ellipsoid 

[64] 


1 

thin plus thick 

mr 


2 

thin plus thick 

m 

disc 

3 

thin plus thick plus halo 

m 


4 

thin plus thick plus halo 



5* 

single maximal disc 

m 


1* 

Ha, HI, HII 

"W 

gas 

2 

Ha, HI, HII 

m\ 


TABLE I: Summary of all models of stellar bulge, stellar 
disc and gas used to describe the baryonic component of our 
Galaxy. For further details, please see Ref. m- The configu¬ 
rations defining the representative baryonic model used later 
on in our analysis are indicated with an asterisk. 


three-dimensional morphologies of each baryonic compo¬ 
nent, namely stellar bulge, stellar disc(s) and gas. Each 
component is now discussed in turn. 

Stellar bulge. The bulge dominates the inner 2 — 
3 kpc of our Galaxy and presents a triaxial shape with 
a bar extending at positive Galactic longitudes. This 
general picture is consistently painted by different obser¬ 
vations, but the morphological details are rather uncer¬ 
tain. We implement seven alternative data-based con¬ 
figurations for the bulge: exponential E2 [59], gaussian 
G2 [55] , gaussian plus nucleus [SD] , truncated power laws 
ElliS], power law plus long bar |62j and double ellip¬ 
soid [63] . All models are normalised to the observed mi- 
crolensing optical depth towards {£, b) = (1.50°, —2.68°), 
(r) = 2.171°;"^ X 10-6 [35]El]. 

Stellar disc(s). The structure of the Galactic disc 
has been thoroughly studied using different photometric 
surveys. The disc is usually modelled as a flattened com¬ 
ponent with a fastly decaying profile and is often sub¬ 
divided into thin and thick populations. As with the 
bulge, there is no consensus over the global morphologi¬ 
cal details. We take therefore hve alternative morpholo¬ 
gies: pure thin plus thick discs [55] [55], thin plus thick 
discs with a stellar halo component [57] [55] and a sin¬ 
gle maximal disc [21]. All models are normalised to the 
recent measurement of the local total stellar surface den¬ 
sity, E* = 38 ± 4 Mq /pc^ [33] . 

Gas. The gas component takes the form of molecular, 
atomic and ionised hydrogen (as well as a small frac¬ 
tion of heavier elements), and its distribution is fairly 
well-known. We model the gas in the inner 10 pc (ef¬ 
fectively, a point-like distribution for our purposes) ac¬ 
cording to Ref. [73] and in the inner 2 kpc according 
to Ref. [73]. Above 2 kpc we implement two alterna¬ 
tive morphologies based on Refs. [551 EO] . The un¬ 
certainties associated to each gas model stem mainly 
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FIG. 1: The rotation curve of our Galaxy. The compilation of rotation curve measurements discussed in the text is shown by 
the red data points, whereas the bracketing of the 70 baryonic models implemented is spanned by the grey band. Both elements 
are plotted with Icr uncertainties. The mean rotation curve predicted by a representative baryonic model (consisting of bulge 
1 [59] . disc 5 [21] and gas 1 [69], cf. Tab. is denoted by the black solid line. For convenience, the angular circular velocity 
is shown in the upper panel and the actual circular velocity in the lower panel. Here we take Ro = 8 kpc, vo = 230km/s and 
(U, V, W)q = (11.10,12.24, 7.25) km/s |57|. 


from the poorly constrained CO-to-H 2 conversion fac¬ 
tor: (2.5— 10) X 10^®cm“^(K km/s)~^ for R < 2kpc and 
(0.5 —3.0) X 10^° cm“^(K km/s)~^ for R > 2kpc [5^171] . 

Tab. H] summarises all models specified above. For 
each bulge, disc and gas model, it is now possible to 
compute the corresponding contribution to the rotation 


curve. With seven bulges, five discs and two gas con¬ 
figurations, there are in total 70 different baryonic dis¬ 
tributions, which we use to bracket the overall baryonic 
contribution -f -b 

In Fig. we show the rotation curve measurements and 
the expected baryonic contribution. The red data points 
show the compilation of kinematic tracers with Icr uncer- 
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FIG. 2: Breakdown of the relative contribution to the overall 
uncertainty on the residual as a. function of Galacto- 

centric radius R. From top to bottom, the bars show the mean 
contribution due to rotation curve data, bulge, disc, gas and 
baryonic bracketing. The larger radial bins at i? > 15 kpc 
are chosen due to the fewer Wc measurements in that re¬ 
gion. As in Fig. we take Rq = 8 kpc, vq = 230km/s and 
([/, V, W)q = (11.10,12.24, 7.25) km/s [S^. 


tainties. The grey band represents instead the envelope 
of the 70 baryonic models, including their corresponding 
la uncertainties, and the black solid line indicates the 
mean curve of a representative baryonic model defined 
by bulge 1, disc 5 and gas 1 (i.e. E2 bulge |59], single 
maximal disc [H] and gas from Ref. [5S], cf. Tab.[l]). The 
upper and lower panels show, respectively, the angular 
circular velocity ujc, which will be used throughout our 
analysis, and the actual circular velocity Vc, which is the 
familiar quantity commonly displayed in rotation curve 
studies. Let us notice once again that the errors of R and 
Vc are strongly (positively) correlated, while those of R 
and ujc are uncorrelated. It is apparent from Fig. [l] that 
baryons fall short of supporting the rotation curve in the 
Milky Way; a detailed analysis treating carefully statis¬ 
tical and systematic uncertainties shows the discrepancy 
is significant already inside the solar circle [53]. 

There are five sources of uncertainty that affect the 
comparison between Wc and and the measurement of 
ojc itself, the normalisations of bulge, disc and gas and 
the baryonic model bracketing. It is instructive to quan¬ 
tify the relative importance of each uncertainty on, for 
instance, the residual — Wb- This is shown in Fig. 
We find that for R < Rq baryonic bracketing and, to 
a lesser extent, the bulge normalisation are the main 
sources of uncertainty, whereas for R > Rq the contri¬ 
bution of rotation curve measurements dominates over 
all other components. The uncertainties on disc and gas 
normalisations are small across the full radial range. In 
order to improve the dynamical constraints on the mass 
distribution in the Milky Way (including the ones shown 
later on in this work), it is therefore crucial to improve 
our knowledge of the spatial distribution of baryons at 
R Ro and to increase the precision of rotation curve 


data at R > Rq. Fortunately, Gaia |75|, APOGFF-2 
(SDSS-IV) [76], WFIRST [77], WFAVF [78] and 4MOST 
[75] should be able to help out on both fronts over the 
coming years. 


III. METHODOLOGY 


Having introduced the observed rotation curve Wc and 
the baryonic contribution oj-q in Sec. |ll| we now proceed 
to interpret their discrepancy in terms of dark matter. 
The goal here is to determine the best-fit parameters for 
a dark matter distribution by fitting the resulting total 
expected curve (baryons and dark matter) to the obser¬ 
vational data. (For a different procedure to extract the 
dark matter profile directly from the data, see Ref. [80]). 
We make the common assumption of a spherical dark 
matter distribution, for which the additional contribu¬ 
tion to the rotation curve simply traces the total dark 
mass enclosed within a given radius R, namely 


2 

^dm 


GMQcn{<R) 

R3 





dr47rr^Pdm(r). 


( 2 ) 


This then leads to a total angular velocity ujt = 
to be compared to the observed uJc- For 
the (spherical) dark matter density in Fq. ([^, we use ei¬ 
ther the generalised Navarro-Frenk-White (NFW) profile 
[5THM]. Pdm oc (r/rs)“^(l -I- r/rs)“^+'^, where 7 is the 
inner slope and Vs is the scale radius, or the Finasto pro¬ 
file [SlllHS], pdm oc exp(— 2 ((r/rs)“ — l)/a), where a is a 
shape parameter. In both cases, we parametrise the nor¬ 
malisation of the profiles in terms of the the local dark 
matter density po = PdmiRo)- 

The remaining task is to fit u>t to oJc and derive con¬ 
straints on the dark matter parameters po, 7 (or a) and 
Ts ■ A subtlety here concerns the large errors on R present 
in the compilation of kinematic data (cf. Fig. [^, which 
preclude the use of a standard chi-square fitting. In¬ 
stead of binning the data in R (and thus necessarily 
losing information), we follow Ref. [53] to define a chi- 
square that accounts for uncertainties on both the verti¬ 
cal and horizontal axes. Defining the dimensionless vari¬ 
ables {x,y) = [R/Rq, Wc/wq — 1) with ujq = vq/Rq, we 
have 


N N 





yt,i) 


2 


-I- a 


2 

h,i 





(3) 


where (xi ± ax,i,yi ± <Xy^i) are the rotation curve mea¬ 
surements discussed in Sec. |llj (Tb,i is the uncertainty 
of the individual baryonic model evaluated at Xi and 
{xt^i, yt,i) are the points that minimise di along the curve 
yt{x) = uJt{R = xRq)/ujq — 1. The sum runs over all 
the N objects in the compilation at i? > i?cut = 2.5 kpc 
[N = 2687 with 2081, 506 and 100 from gas, stars and 
masers, respectively) in order to exclude the innermost 
regions of the Galaxy where axisymmetry breaks down 
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FIG. 3: Dark matter fit to the observed rotation curve of our Galaxy. The left (right) panels show the favoured regions of 
the dark matter parameter space for a generalised NFW (Einasto) profile with a fixed scale radius rs = 20kpc. The upper 
panels display with grey lines the 2cr confidence regions corresponding to each baryonic model. For the representative baryonic 
model | 21l 1591169 | , besides the 2a confidence region in thick black, we also plot the 5a goodness-of-fit region in thin black 
dashed. The bottom panels show for each baryonic model the profiled local dark matter density range encompassed by the 2a 
confidence regions (thick) and by the 5a goodness-of-fit regions (thin). The baryonic models are ordered from top to bottom 
as in Tab. [n| The representative model is indicated by the arrows. This figure assumes Ro = 8kpc, vo = 230km/s and 
(U,V,W)q = (11.10,12.24, 7.25) km/s [HZ|. For reference, 0.38GeV/cm® = 0.01Mq/pc®. 


and some tracers may present non-circular orbits. We 
note in passing that Eq. ([^ can be applied since x and 
y (i.e. R and ujc) have uncorrelated errors. All results 
shown in the remainder of the work are based on the chi- 
square statistic defined above. We wish to point out that 
in Ref. [33] we have checked through Monte Carlo cal¬ 
culations that this statistic follows closely a chi-square 
distribution with the appropriate number of degrees of 
freedom for the case of our representative baryonic model 
and using the typical radial and velocity uncertainties 
found in our kinematic compilation. 


IV. RESULTS 

We now turn to deriving constraints on the Galactic 
dark matter profile. For each of the 70 baryonic models 
described in Sec. |ll| a scan was performed over the dark 
matter profile parameter space and the computed as 


in Eq. We first focus our attention on the parame¬ 
ter space {po, 7 ) for NFW and {po, a) for Einasto ignoring 
the correlation between these parameters and other quan¬ 
tities such as Tg, Rq, Vq or Fq, and then we explore quan¬ 
titatively the effect of varying the scale radius. Galactic 
parameters, data selection and systematics. 

Fig.i presents the main results of our analysis for 
the generalised NFW (left) and Einasto (right) profiles, 
both with scale radius = 20 kpc and for our baseline 
Galactic parameter configuration with i?o = 8 kpc, Vq = 
230 km/s and(C/,U,W)Q = (11.10,12.24, 7.25) km/s [S7]. 
We convey our results with two distinct regions in the 
parameter space: (i) the 95.45% (2cr) confidence region 
encompassed by < Xbf + where Xbf is the x^ 
of the best fit configuration and = 6.18 (corre¬ 

sponding to two fitted parameters at 2a); and (ii) the 
5a goodness-of-fit region defined by < xicry where xia 
is the y^ corresponding to an equivalent 5cr significance 
(i.e. p-value 2.87 x 10“^). The top panels display the 2cr 
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FIG. 4: The dependence of the dark matter ht on the scale 
radius. The plot shows the favoured regions of the dark mat¬ 
ter parameter space for a pure NFW prohle (7=1) and 
the representative baryonic model |21[ 1591 I69j . The con¬ 
tours and line coding are the same as in Fig. (top left). 
We have assumed here Ro = 8kpc, vo = 230km/s and 
(U,V,W)q = (11.10,12.24, 7.25) km/s [57]. For reference, 
0.38GeV/cm^ = O.OIMq/pc^. 

confidence regions for all baryonic models and also the 
5(7 goodness-of-fit region for the representative baryonic 
model [211 IMl ES] ■ The bottom panels show instead the 
profiled range of local dark matter density encompassed 
by the 2cr confidence region and the 5a goodness-of-fit 
region of all baryonic models. In Tab. [ll] we report the 
best fits for all baryonic models along with the profiled 
ranges of local dark matter density. 

The constraints in Fig. are significantly more sensi¬ 
tive to the local dark matter density than to the slope 
of the profile. This simply reflects the impact of these 
two quantities on the dark matter content of the Galaxy 
and thus on its contribution to the rotation curve. Over¬ 
all, the remarkably precise 2a confidence regions arise 
from the choice of an unbinned analysis and the large 
number of kinematic tracers adopted. However, Fig. 
clearly shows how confidence regions corresponding to 
different baryonic models shift sizeably across the param¬ 
eter space, a feature which reflects the importance of the 
baryonic contribution to the rotation curve of the inner 
Galaxy. In other words, the precision currently allowed 
by kinematic measurements is hindered by large system- 
atics associated with the distribution of baryons. This 
highlights the reason for adopting a comprehensive collec¬ 
tion of baryonic models in our analysis: whereas a single, 
fiducial model may well represent a useful benchmark, it 
cannot be used for accurate dynamical constraints unless 
one has an a priori expectation about the actual baryon 
morphology of the Galaxy. The spread of the dark mat¬ 
ter constraints over the whole range of baryonic models is 
a measure of the systematic uncertainty due to baryonic 
modelling, which cannot be dealt with by averaging nor 
by marginalising over the entire range of allowed mor¬ 
phologies. We therefore present the constraints for all 


baryonic models instead of focussing on a single model 
with unknown systematics. 

Another interesting message conveyed by Fig. is that 
dynamical constraints start to be able to resolve be¬ 
tween different baryonic models. Let us focus on the 
left panels corresponding to the NFW profile (very simi¬ 
lar considerations can be drawn for the Einasto profile). 
From the bottom left plot it is clear that there are five 
blocks of baryonic models clustered around slightly dif¬ 
ferent local dark matter densities. These five blocks cor¬ 
respond precisely to the five discs implemented in our 
analysis (cf. Tab.|l]). The inferred po ranges from around 
0.50 GeV /eve? for models including disc 1 down to around 
0.42 GeV /env for models with disc 5, and in between for 
the other discs 2-4. These figures are comparable to but 
somewhat higher than values typically found when ap¬ 
plying global methods [23 dS] ■ The inner slope is less 
well determined, but the confidence regions in the top 
left panel of Fig. are mainly centred around three val¬ 
ues of the inner slope (7 ^ 0.6,1.2,1.4) with a few out¬ 
liers. This configuration correlates with the bulge used 
in each baryonic model. When using any of the discs 1- 
4, the baryonic models with the bulges 1-4 or 6 point to 
7 ~ 1.2, the ones with bulge 5 cluster around 7 ^ 0.6 
and the ones with bulge 7 tend towards 7 ~ 1.4. The 
case of models with disc 5 is slightly different: the bulges 
1-4 or 6 point now to 7 ~ 0.6, while bulge 5 prefers shal¬ 
low profiles 7 ^ 0 (in this case cored profiles are possibly 
favoured over NFW or Einasto) and bulge 7 points to 
7 ^ 1 . From all this discussion, we learn that the lo¬ 
cal dark matter density is particularly sensitive to the 
morphology of the stellar disc, whereas the inner slope of 
the dark matter profile is mostly dependent on the bulge 
configuration. The gas plays a relatively minor role. 

It is worth noticing that there are eight baryonic mod¬ 
els for which the 5cr goodness-of-fit region is vanishing 
(cf. bottom panels in Fig. and Tab. 0. This is valid 
for both the generalised NFW and the Einasto profiles. 
These eight baryonic models coupled to all the dark mat¬ 
ter profile configurations tested fall short of supporting 
the observed rotation curve and are excluded at 5a level 
or more. The common characteristic of the eight mod¬ 
els is the use of bulge 7 coupled to any of the two gas 
configurations and any of the discs except disc 5. If one 
were to have independent evidence for the dark matter 
profiles tested in Fig. then bulge 7 (i.e. the double el¬ 
lipsoid bulge [64]) would be too light and disfavoured by 
current dynamical data unless it is associated with disc 
5 (i.e. the single maximal disc [H]). 

For the representative baryonic model EH [Ml ES 
(highlighted in black in the panels of Fig. [^, the 2a 
ranges for the local dark matter density read 

NFW: Po = 0.420j:^;)Ji8 {2a) ± 0.025 GeV/cm^ , (4) 

Einasto: po = 0.420;)! )J 2 ? (2cr) ± 0.026 GeV/cm^ , (5) 

where the second error is the standard deviation of the 
best-fit values over all baryonic models, which we use here 
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FIG. 5: The dependence of the dark matter fit on Galactic parameters, data selection and systematics. In the left panel 
we show the effect of different values of Ro, vq and Vq, of a systematic due to spiral arms and of an enhanced stellar disc 
normalisation uncertainty. The right panel illustrates the impact of a different data selection (gas, stars, masers) as well as the 
result of a standard binned analysis for comparison. In all cases the 2a confidence region for the generalised NFW profile and 
the representative baryonic model is plotted. The thick black contour is the baseline case presented in Fig. 0 (top 

left). For reference, 0.38GeV/cm^ = O.OIMq/pc^. 


as a measure of the current systematic uncertainty due 
to baryonic modelling. The interested reader can find in 
Tab. |TT] the inferred ranges for all baryonic models. For 
the sake of completeness, let us also point out that all val¬ 
ues reported assume a spherical dark matter profile; an 
oblate profile as suggested by numerous numerical simu¬ 
lations would lead to a higher local dark matter density 
by tens of percent [HHl HZ] . 

We now comment on the impact of breaking the as¬ 
sumptions behind our main results shown in Fig. The 
effect of letting free the scale radius is illustrated in Fig.j^ 
for a pure NFW profile with inner slope 7 = 1 and the 
representative baryonic model HU US HD- For compar¬ 
ison, numerical simulations of Milky Way-like halos find 
scale radii in the range ~ 12 — 37kpc [55H5D] . As the 
kinematic data used in the current paper are restricted to 
the inner Galaxy (namely, i? < 20kpc), the 5 ct exclusion 
region is rather loose unless Vg < 20kpc. Better con¬ 
straints on the scale radius require tracers in the outer 
Galaxy, as pursued e.g. in Ref. HHIHI]. Notwithstanding, 
Fig. HI clearly shows that the local dark matter density 
is well constrained with uncertainties comparable to the 
ones obtained in Fig. (left) where we fixed Vg = 20 kpc 
but varied 7. 

Finally, Fig. [^depicts the impact of Galactic parame¬ 
ters, data selection and systematics in our dynamical con¬ 
straints. For this figure we only plot the 2a confidence re¬ 
gion for the generalised NFW profile and the representa¬ 
tive baryonic model HU EH] ES] . The left panel shows the 
effect of changing the Galactic fundamental parameters 
Rq, Vq and Vq. The current ranges for these quantities 
read Rq = 8.0 ± 0.5kpc [5^ ICT155] . vq = 230 ± 20km/s 
^ |94H97| and Fq = 5.25 - 26km/s [5l|53 |92||9H|, 
but there are important correlations among the differ¬ 
ent values. Let us consider specifically the measurement 
Rq = 8.33 ± 0.35 kpc HI] from the monitoring of stellar 


orbits around the central supermassive black hole and 
Qq = = 30.26 ± 0.12km/s/kpc |99| from the 

proper motion of Sgr A*. Put together, the two mea¬ 
surements define four Icr configurations depending on the 
adopted value of Vq : 

(a) Rq = 7.98 kpc, Vq = 214.52 km/s for Vq = 26km/s; 

(b) Ro = 7.98 kpc, vq = 237.18 km/s for Vq = 5.25 km/s; 

(c) Rq = 8.68 kpc, Vo = 235.62 km/s for Vq = 26 km/s; 

(d) Ro = 8.68 kpc, vg = 258.45 km/s for Vq = 5.25 km/s. 

In the left panel of Fig. we present the results for the 
two most extreme configurations, which turn out to be 
(a) and (d). Glearly, the current (correlated) uncertain¬ 
ties on Ro, Vq and Vq hinder the determination of the 
local dark matter density, while their effect on the deter¬ 
mination of the inner slope is somewhat less important. 
Also shown in the left panel is the impact of the sys¬ 
tematic motion due to spiral arms modelled according 
to Ref. [H] and of an enhanced 20% uncertainty on the 
local total stellar surface density (see the supplementary 
information of Ref. [33] for further details). Our con¬ 
straints are fairly robust in both cases. The right panel 
of Fig. id isplays instead how our confidence regions shift 
selecting separately gas kinematics, star kinematics and 
masers. Although there is some dependence on data se¬ 
lection, we opt to show our main results including all data 
available in the compilation of kinematic measurements. 
For comparison, we also plot the results of a standard 
binned analysis applied to the full data set, which is less 
precise overall but more robust against data selection. 


V. CONCLUSION 

In the cold dark matter paradigm, the innermost re¬ 
gions of the Milky Way are expected to harbour a signif- 





icant amount of dark matter. Testing such expectation 
with the help of observations has historically been diffi¬ 
cult due to the uncertainties on both the distribution of 
baryons and the rotation curve. We have tried here to 
overcome those difficulties by combining a comprehensive 
compilation of rotation curve measurements with state- 
of-the-art baryonic models. This allowed us to effectively 
subtract the contribution of baryons off the observed ro¬ 
tation curve, and to identify the favoured range of the pa¬ 
rameters of different dark matter profiles. We believe the 
data and analysis techniques presented here will be use¬ 
ful for future dynamical studies extending beyond stan¬ 
dard assumptions (e.g. non-spherical dark matter pro¬ 
files). At present, the precision allowed by kinematic 
measurements is entirely overshadowed by the current 
uncertainty on the fundamental Galactic parameters and 


on the morphology of baryons. In that respect, the forth¬ 
coming data from the Gaia mission and surveys such as 
APOGEE-2 (SDSS-IV), WFIRST, WEAVE and 4MOST 
will play a crucial role in shrinking these uncertainties 
and will hopefully open a precision era in the measure¬ 
ment of the dark matter distribution in the Galaxy. 
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generalised INl^'W = 2Ukpcj 

Einasto = 2Ukpc) 

baryonic model 

Xhf/^ PO [GeV/cm^] 

Xbf/^ PO [GeV/cm^] 

1 

1-1-1 


1.016 0.479-0.509 (0.386-0.614) 

1.012 0.479-0.510 (0.392-0.614) 

2 

1-1-2 


1.042 0.483-0.514 (o.402-0.604) 

1.037 0.484-0.514 (0.404-0.604) 

3 

2-1-1 


1.044 0.479-0.511 (0.401-0.599) 

1.040 0.482-0.510 (o.402-0.599) 

4 

2-1-2 


1.071 0.484-0.513 (o.418-0.585) 

1.067 0.485-0.515 (o.422-0.585) 

5 

3-1-1 


0.899 0.479-0.514 (o.329-0.686) 

0.893 0.479-0.514 (o.336-0.688) 

6 

3-1-2 


0.919 0.482-0.517 (o.339-0.680) 

0.913 0.483-0.516 (o.345-0.681) 

7 

4-1-1 


1.020 0.480-0.511 (0.387-0.614) 

1.015 0.481-0.511 (0.391-0.615) 

8 

4-1-2 


1.046 0.484-0.514 (o.405-0.604) 

1.041 0.484-0.514 (o.412-0.604) 

9 

5-1-1 


0.744 0.479-0.522 (o.240-0.730) 

0.737 0.478-0.521 (o.245-0.750) 

10 

5-1-2 


0.757 0.482-0.524 (o.248-0.734) 

0.751 0.480-0.523 (o.252-0.753) 

11 

6-1-1 


0.940 0.478-0.512 (o.347-0.661) 

0.934 0.479-0.511 (o.351-0.662) 

12 

6-1-2 


0.962 0.482-0.515 (o.360-0.653) 

0.956 0.482-0.515 (o.366-0.655) 

13 

7-1-1 


1.359 0.476-0.499 

1.361 0.482-0.502 

14 

7-1-2 


1.405 0.482-0.504 

1.406 0.488-0.507 

15 

1-2-1 


1.031 0.453-0.483 (0.371-0.578) 

1.028 0.456-0.483 (0.378-0.575) 

16 

1-2-2 


1.057 0.457-0.487 (o.387-0.565) 

1.054 0.461-0.488 (o.391-0.566) 

17 

2-2-1 


1.061 0.454-0.483 (o.387-0.559) 

1.059 0.456-0.484 (o.390-0.558) 

18 

2-2-2 


1.089 0.458-0.487 (o.407-0.543) 

1.087 0.461-0.488 (o.413-0.545) 

19 

3-2-1 


0.914 0.455-0.490 (o.315-0.652) 

0.910 0.456-0.489 (o.326-0.653) 

20 

3-2-2 


0.936 0.459-0.492 (o.327-0.645) 

0.931 0.460-0.492 (o.334-0.646) 

21 

4-2-1 


1.034 0.453-0.484 (o.370-0.577) 

1.031 0.456-0.484 (o.376-0.576) 

22 

4-2-2 


1.061 0.458-0.488 (o.390-0.564) 

1.057 0.460-0.488 (o.398-0.565) 

23 

5-2-1 


0.762 0.458-0.500 (o.229-0.702) 

0.756 0.457-0.499 (o.235-0.720) 

24 

5-2-2 


0.776 0.461-0.502 (o.238-0.705) 

0.770 0.459-0.501 (o.244-0.723) 

25 

6-2-1 


0.956 0.454-0.486 (o.332-0.626) 

0.952 0.455-0.487 (o.339-0.625) 

26 

6-2-2 


0.979 0.458-0.490 (o.347-0.617) 

0.975 0.460-0.490 (o.351-0.618) 

27 

7-2-1 


1.382 0.447-0.468 

1.386 0.453-0.469 

28 

7-2-2 


1.429 0.454-0.474 

1.434 0.459-0.475 

29 

1-3-1 


1.031 0.452-0.481 (0.370-0.576) 

1.029 0.454-0.481 (0.377-0.573) 

30 

1-3-2 


1.058 0.455-0.486 (o.385-0.563) 

1.055 0.459-0.486 (o.390-0.564) 

31 

2-3-1 


1.062 0.453-0.482 (o.385-0.557) 

1.060 0.454-0.484 (o.389-0.557) 

32 

2-3-2 


1.090 0.457-0.486 (0.406-0.540) 

1.087 0.459-0.486 (o.411-0.542) 

33 

3-3-1 


0.915 0.454-0.488 (o.313-0.650) 

0.910 0.455-0.488 (o.325-0.651) 

34 

3-3-2 


0.936 0.457-0.491 (o.326-0.643) 

0.931 0.458-0.491 (o.333-0.644) 

35 

4-3-1 


1.035 0.452-0.482 (o.370-0.576) 

1.031 0.454-0.483 (o.374-0.574) 

36 

4-3-2 


1.061 0.456-0.486 (o.388-0.562) 

1.058 0.459-0.486 (o.397-0.563) 

37 

5-3-1 


0.762 0.456-0.498 (o.228-0.701) 

0.756 0.455-0.497 (o.234-0.719) 

38 

5-3-2 


0.777 0.459-0.501 (o.237-0.704) 

0.770 0.458-0.500 (o.242-0.722) 

39 

6-3-1 


0.956 0.452-0.485 (o.331-0.624) 

0.952 0.453-0.485 (o.338-0.623) 

40 

6-3-2 


0.979 0.456-0.488 (o.345-0.615) 

0.975 0.458-0.489 (o.350-0.616) 

41 

7-3-1 


1.383 0.445-0.466 

1.387 0.453-0.469 

42 

7-3-2 


1.430 0.452-0.472 

1.435 0.459-0.475 

43 

1-4-1 


0.960 0.428-0.458 (0.317-0.588) 

0.959 0.431-0.459 (0.329-0.586) 

44 

1-4-2 


0.982 0.432-0.462 (o.327-0.582) 

0.981 0.434-0.464 (o.343-0.579) 

45 

2-4-1 


0.985 0.429-0.459 (o.326-0.576) 

0.984 0.431-0.459 (o.342-0.573) 

46 

2-4-2 


1.008 0.433-0.462 (o.341-0.567) 

1.007 0.435-0.463 (o.349-0.567) 

47 

3-4-1 


0.844 0.429-0.464 (o.266-0.660) 

0.842 0.430-0.464 (o.281-0.661) 

48 

3-4-2 


0.861 0.433-0.467 (o.275-0.655) 

0.859 0.434-0.468 (o.290-0.655) 

49 

4-4-1 


0.966 0.430-0.461 (o.319-0.590) 

0.964 0.431-0.462 (o.325-0.589) 

50 

4-4-2 


0.988 0.434-0.465 (o.331-0.583) 

0.986 0.435-0.466 (o.338-0.581) 

51 

5-4-1 


0.689 0.437-0.481 (o.186-0.717) 

0.684 0.436-0.480 (o.201-0.735) 

52 

5-4-2 


0.701 0.440-0.483 (o.193-0.721) 

0.696 0.439-0.482 (o.209-0.738) 

53 

6-4-1 


0.884 0.428-0.462 (o.284-0.635) 

0.882 0.431-0.462 (o.295-0.634) 

54 

6-4-2 


0.902 0.432-0.465 (o.292-0.629) 

0.901 0.434-0.465 (0.301-0.629) 

55 

7-4-1 


1.289 0.422-0.444 

1.293 0.429-0.446 

56 

7-4-2 


1.327 0.428-0.449 

1.331 0.435-0.452 

57* 

1-5-1 


0.773 0.402-0.441 (0.199-0.656) 

0.772 0.399-0.439 (0.208-0.671) 

58 

1-5-2 


0.783 0.405-0.444 (0.207-0.659) 

0.782 0.402-0.441 (o.217-0.674) 

59 

2-5-1 


0.785 0.402-0.441 (o.205-0.656) 

0.784 0.399-0.438 (o.215-0.669) 

60 

2-5-2 


0.796 0.404-0.443 (o.213-0.658) 

0.795 0.403-0.441 (o.221-0.670) 

61 

3-5-1 


0.690 0.403-0.446 (o.161-0.672) 

0.689 0.401-0.445 (0.183-0.689) 

62 

3-5-2 


0.699 0.406-0.448 (o.168-0.676) 

0.698 0.404-0.446 (0.189-0.694) 

63 

4-5-1 


0.784 0.411-0.451 (0.207-0.633) 

0.782 0.409-0.449 (o.212-0.651) 

64 

4-5-2 


0.795 0.414-0.453 (o.215-0.637) 

0.793 0.412-0.451 (o.220-0.655) 

65 

5-5-1 


0.576 0.411-0.441 (o.098-0.619) 

0.576 0.411-0.455 (0.117-0.641) 

66 

5-5-2 


0.583 0.415-0.449 (o.102-0.626) 

0.582 0.413-0.461 (0.121-0.648) 

67 

6-5-1 


0.714 0.403-0.445 (o.173-0.668) 

0.713 0.401-0.443 (0.191-0.685) 

68 

6-5-2 


0.724 0.405-0.447 (0.179-0.672) 

0.723 0.403-0.445 (0.195-0.688) 

69 

7-5-1 


1.007 0.389-0.419 (o.295-0.529) 

1.008 0.385-0.416 (o.296-0.525) 

70 

7-5-2 


1.022 0.393-0.424 (0.305-0.526) 

1.022 0.391-0.420 (0.309-0.520) 


TABLE II: The best fit and local dark matter density inferred using generalised NEW and Einasto profiles with fixed scale 
radius Vs — 20kpc for all baryonic models. The baryonic models are specified by the configurations of bulge, disc and gas 
(cf. Tab.[^ in the form [b-d-g], where b,d,g stand for bulge, disc and gas, respectively. Besides the best fit for each model, 
we also report the profiled ranges of local dark matter density covered by the 2a conhdence region (no parentheses) and by 
the 5 ct goodness-of-fit region (between parentheses) shown in the bottom panels of Fig. The representative baryonic model 
is indicated with an asterisk. Here we take Ro = 8kpc, uo = 230km/s and {U,V,W)q — (11.10,12.24, 7.25) km/s |57| . For 
reference, 0.38GeV/cm^ = O.OIMq/pc^. 
















